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We investigate in detail the momentum dependence {MD) of the ef- 
fective in medium Nucleon-Nucleon (NN) interaction in the isovec- 
tor channel. We focus the discussion on transport properties of 
the expected neutron-proton (n/p) effective mass splitting at high 
isospin density. We look at observable effects from collective flows 
in Heavy Ion Collisions {HIC) of charge asymmetric nuclei at in- 
termediate energies. Using microscopic kinetic equation simulations 
nucleon transverse and elliptic collective flows in Au + Au collisions 
are evaluated. In spite of the reduced charge asymmetry of the in- 
teracting system interesting isospin — MD effects are revealed. 
Good observables, particularly sensitive to the n/p-mass splitting, 
appear to be the differences between neutron and proton flows. 
The importance of more exclusive measurements, with a selection 
of different bins of the transverse momenta [pt) of the emitted par- 
ticles, is stressed. In more inclusive data a compensation can be 
expected from different pt-contributions, due to the microscopic 
isospin — MD structure of the nuclear mean field in asymmetric 
matter. 
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1 Introduction 

A key question in the physics of unstable nuclei is the knowledge of the EOS 
for asymmetric nuclear matter away from normal conditions. We remind again 
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the effect of the symmetry term at low densities on the neutron skin struc- 
ture, while the knowledge in high densities region is crucial for supernovae 
dynamics and neutron star cooling. The paradox is that while we are plan- 
ning second and third generation facilities for radioactive beams our basic 
knowledge of the symmetry term of the nuclear Equation of State EOS is 
still extremely poor. Effective interactions are obviously tuned to symmetry 
properties around normal conditions and any extrapolation can be quite dan- 
gerous. Microscopic approaches based on realistic A^A^ interactions, Brueckner 
or variational schemes, or on effective field theories show a quite large variety 
of predictions, see the discussions in [1,2] and refs. therein. 

In this paper we look at the mean field momentum dependence (MD) in 
the isospin channel and study its dynamical effects in heavy ion collisions 
at intermediate energies. This represents a completely open problem as we 
can quickly realize from a simple analysis of the effective Skyrme-like non- 
relativistic forces. Let's look at the local and non-local contributions to the 
potential symmetry term of various Skyrme interactions, SHI, SGII, SKM* 
(see [3] and refs. therein) and the more recent Skyrme-Lyon forms, SLy230a 
and Sly230b (or SlyA), [4]. We can clearly see a sharp transition from the 
previous Skyrme forces to the Lyon parametrizations, with an almost inver- 
sion of the signs of the two contributions. The important repulsive non-local 
part of the Lyon forces leads to a completely different behavior of the neutron 
matter EOS, of great relevance for the neutron star properties. Actually this 
substantially modified parametrization was mainly motivated by a very un- 
pleasant feature in the spin channel of the previous Skyrme forces, the collapse 
of polarized neutron matter, see discussion in [5,4]. 

A related interesting effect can be seen on the neutron/proton effective masses. 
For an asymmetry P = ^^-^ = 0.2 (the ^^'^Au asymmetry) we can get a 
splitting up to the order of 15% at normal density po, and increasing with 
baryon density. 

The sign itself of the splitting is very instructive. Passing from Skyrme to 
Skyrme-Lyon we see a dramatic inversion in the sign of the n/p mass split- 
ting. The Lyon forces predict in n-rich systems a neutron effective mass al- 
ways smaller than the proton one. We will come back to this point. Here we 
just note that the same result is obtained from more microscopic relativistic 
Dirac-Brueckner calculations [6] and in general from the introduction of scalar 
isovector virtual mesons in Relativistic Mean Field [RMF) approaches [7,8]. 
At variance, non-relativistic Brueckner-Hartree-Fock calculations are leading 
to opposite conclusions [9,10]. This will be one of the main questions to address 
in the reaction dynamics of exotic nuclear systems. 

We can expect important effects on transport properties ( fast particle emis- 
sion, collective flows, resonance and particle production around the threshold) 
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of the dense and asymmetric NM that will be reached in Radioactive Beam 
collisions at intermediate energies. 

This is the aim of the present paper. We focus our attention on the isospin 
dependence of collective flows. We perform collision simulations based on mi- 
croscopic kinetic equation models. We start from realistic effective interactions 
widely used for symmetric systems. We test very different parametrizations 
in the momentum dependence of the isovector channel, taking care that the 
symmetry energy, including its density dependence, will be not modified. We 
study in particular the transport effect of the sign of the n/p effective mass 
splitting m* — m* in asymmetric matter at high baryon and isospin density. 

We see very interesting isospin — MD effects on n/p transverse and elliptic 
flows in semi central Au + Au collisions at 250 AMeV beam energy. The impor- 
tance of more exclusive data, with a good selection of the transverse momenta 
of the emitted particles, is stressed. 

In Sect. 2 we present details of the chosen effective interactions. In Sect. 3 the 
collective flow results are discussed and finally in Sect. 4 some conclusions are 
drawn with relative perspectives. 



2 The choice of the effective interactions 



The phenomenological non-relativistic effective EOS used in this work was 
first introduced by I.Bombaci et al. [11,12] in a general form suitable for as- 
trophysical and heavy ion physics applications. We have modified the isovector 
part in order to get two new parametrizations with different n/p effective mass 
splittings while keeping the same symmetry energy at saturation, including a 
very similar overall density dependence. 

We start writing down the energy density as a function of the asymmetry 
parameter /3 = "^^^^i 
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£MD{Pn, Pp, P) = C — {In + Xp) + ^— /3(J„ - Ip) 

Po 5 p„ 

where the integrals Xr{k,A) = J d^kfr{k)g{k, A) include the momentum 

r 21 ~i 

dependent part of the mean field g{k,A) = 1 + (^^^) ; the subscript 

T = n,p stands for neutrons and protons respectively; p^ = 0.16 /m~^ is the 
normal density of nuclear matter. We refer to this parametrization as the 
BGBD — EOS. For symmetric nuclear matter {(3 = 0), the energy density 
Eq.(l) reduces to the parametrization proposed by Gale, Bertsch and Das 
Gupta {GBD interaction, [13]). The parameters A, B, C, a and A take the 
same values as in [13] {A = -144 Me V, B = 203.3 Mel/, C = -75MeV, 
0" = |, A = 1.5p^p\ where p^p^ is the Fermi momentum at normal density), 
and provide a soft EOS for symmetric matter, with a compressibility Knm — 
210 MeV. 

The value of zi sets the strength of the momentum dependence (MD) in the 
isospin channel; the remaining parameters Xq and X3 can be set to fix the 
symmetry energy. 

Taking the functional derivative of the energy density with respect to the 
distribution function f^-, we obtain, apart from the kinetic term Ekin, the mean 
field potential for neutrons and protons: 




Fig. 1. Potential symmetry energy as a function of density (solid line) for m* < m* 
(left) and m* > m* (right). Dotted lines refer to density-dependent contributions, 
dashed-dotted to momentum dependent ones. Small panels indicate the correspond- 
ing behaviour of proton and neutron effective masses as a function of density, with 
/3 = 0.2. 
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where the subscripts in the integrals are r ^ r'\ the upper signs refer to 
neutrons, the lower ones to protons. Protons and neutrons in asymmetric 
matter experience different interactions; when /9 = we return to the GBD 
mean field. It is known that a mean field MD implies the idea of an in-medium 
reduction of nucleon mass and leads to the introduction of effective mass [14]. 

The last BGBD term includes, besides the usual GBD momentum depen- 
dence, an isospin-dependent part from which we get different effective masses 
for protons and neutrons. In fact, the effective mass has the following defini- 
tion: 



m 
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and we see that mass splitting is determined not only by different momentum 
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Fig. 2. Mean field potential as a function of momentum at normal density, for an 
asymmetry (3 = 0.2. Solid lines refer to the case m* > m*, dashed lines to m* < m*. 
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dependence of mean field, but also by Fermi momenta of neutrons and protons. 
In our case we get: 



m 



2m 



5 '^J po 



1+ (l±/5)(2/3)(J^V2/3) 



In order to investigate the effects of mass splitting on non-relativistic heavy 
ion collisions we choose two sets of parameters (shown in Table 1) which give 
opposite splitting, but quite similar behaviour of the symmetry energy. 



mass splitting 


Zl 




X3 


m* > m; 


-36.75 


-1.477 


-1.01 


m* < m; 


50 


1.589 


-0.195 



Table 1 

Values of the parameters zi, xq and X3 for opposite mass splitting, giving the same 
Esym{po) = 33AMeV. 

Assuming the usual parabolic law for the EOS of asymmetric nuclear matter: 

we fix the same symmetry energy at the saturation density EgymiPo) ~ 33 AMeV. 

Figure 1 shows the potential part of the symmetry energy as a function of den- 
sity (solid line) for the two choices of mass splitting. We remark the very sim- 
ilar overall density dependence. The dashed-dotted and dotted lines indicate 
respectively the contributions from the momentum dependent and density- 
dependent part of the EOS. A change in the relative sign of mass splitting is 
related to opposite behaviours of these two contributions, exactly as already 
noticed in the Introduction for the Skyrme-like forces. Small panels on the top 
left of each graph illustrate the corresponding mass splitting as a function of 
density for an asymmetry (3 = 0.2 (the ^^"^Au asymmetry). 

Let us discuss now the relation between effective mass and momentum depen- 
dence. From the definition Eq.(3) we see that the effective mass is inversely 
proportional to the slope of mean field at the Fermi momentum. Then, for 
nucleons with smaller effective masses the potential will be more repulsive at 
momenta higher than the Fermi one, and more attractive at low momenta. 
Mean field potential at normal density as a function of momentum is shown 
in Figure 2 for a choice (3 = 0.2. As we see, the parametrizations used in this 
work give rise to opposite behaviors for low and high momentum particles. 
This is a very general feature, due to the meaning itself Eq.(3) of effective 
mass. As we will see in the reaction dynamics this behavior can give rise to 
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some compensation effects between low and high momenta contributions. All 
that will lead to an expected larger isospin — MD sensitivity of more exclu- 
sive measurements, in particular with a tranverse momentum selection of the 
nucleons emitted for a given rapidity. 



3 Effective Masses and Collective Flows 

Collective flows in heavy ion reactions have been widely used to investigate 
properties of nuclear matter, such as compressibility or in-medium cross sec- 
tions [15], isospin dependence [16,17,18] and to find out signals of the mean 
field MD ([19], [20]). The problem of momentum dependence in the isospin 
channel is still open, and recently some possible effects at relativistic ener- 
gies have been studied [21]. Intermediate energies are important in order to 
have high momentum particles and to test regions of high baryon (isoscalar) 
and isospin (isovector) density during the reactions dynamics. Now we present 
some qualitative features of the dynamics in heavy ion collisions related to the 
splitting of nucleon effective masses. 

Our simulations are performed by means of the Boltzmann-Nordheim-Vlasov 
BNV transport code from Ref. [22], implemented with a BGBD — like mean 
field, and using free (isospin dependent) isotropic cross sections in the colli- 
sion term. We focus on semi-peripheral Au + Au collisions at the energy of 
250 AMeV and we evaluate transverse and elliptic flows of neutrons and pro- 
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Fig. 3. Balance energy in Au + Au collisions for intermediate impact parameters. 
Solid circles refer to the case m* > m*, empty circles to m* < m*. Triangles 
represent the data from ref. [25] 
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tons. We consider as emitted particles those resulting in density regions below 
1/8 Pq during the reaction. We have checked that the results are not depending 
on the size of the box where the simulation is performed and on the choice of 
the freeze-out time. 

These flow observables can be seen respectively as the first and second coeffi- 
cients from the Fourier expansion of the azimuthal distribution: 

^{y.Pt) = 1 + Vicos{(P) + 2V2Cos{2<p) 



where pt = ^vl + Vy is the transverse momentum and y the rapidity along 
beam direction. 



3.1 Tranverse flows 



The transverse flow can be expressed as: 

Pt 

It provides information on the azimuthal anisotropy of the transverse nucleon 
emission and has been recently used to study the EOS and cross section 
sensitivity of the balance energy [15]. 

Actually the evaluation of the balance energy for the proton more inclusive 
Directed Flow, integrated over all transverse momenta see [23,16], in Au + Au 
collisions represents a very good test for our effective interactions at lower 
energies. The results are shown in Fig. 3 and compared with NSCL-MSU data. 
Our proton directed flows are calculated for semicentral collisions, i.e. bred = 
0.5 [24]. 

Both parametrizations are well reproducing the experimental balance energy 
region. As expected we see a slightly earlier balance of the flow in the m* > m* 
case due to some more repulsion for fast protons but the effect is hardly 
appreciable at these energies. We will look more carefully at these effects with 
increasing beam energy and for more exclusive observables. 

We will present now n/p transverse flow results for semicentral ^^"^Au -|-^^^ Au 
collisions at 250 AMeV beam energy, for particles in regions of relatively high 
rapidities. It is useful to work with momenta and rapidities normalized to 
projectile ones in the center of mass system (cm), defined as y^^^ = {y/yproj)cm 
and pf'^ = {pt/Pproj)cm- Results of our calculations are illustrated in Fig. 4, 
where the left and right panels illustrate the two cases m* < m* and m* > m* 
respectively. 
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High rapidity selections allow us to find effects of the mean field MD, which 
becomes more important with increasing momenta. Indeed the different be- 
haviour of mean field depending on mass splitting is evident all over the range 
of transverse momentum shown here. We emphasize that the average mo- 
mentum of these particles is generally beyond the Fermi momentum. Indeed a 
rough estimate in the center of mass system gives already for the beam parallel 




"^■^^0 0.5 I 1.5 

Fig. 5. pt dependence of the difference between neutron and proton transverse flows 
in Au + Au collisions, same energy, centrality and rapidity selection as in Fig. 4. 
Solid circles refer to the case m* > m*, empty circles to m* < ni*. 
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component 

p.^mo- 0.8{yproj)cm ^ 260 - 280 MeV/c 
In this momentum region a smaller effective mass determines greater repulsive 
interaction. When m* < m* neutrons feel greater repulsion than protons and 

their deflection in the reaction plane for p[^^ < 0.6 is close to protons, where 
the Coulomb field is also acting. At variance in the case m* > m* the proton 
flow is greatly enhanced if compared with neutrons. 

At high transverse momentum we must take into account the increase of 
the out-of-plane component, so in this case increased repulsion determines 
a greater reduction of transverse flow. It causes a crossing between proton and 
neutron flow in the case m* > m*. This behaviour is clearly seen if we look at 

the different slope of the curves in Fig. 4 (right panel) around pf*^ > 0.6. 

A quite suitable way to observe such effective mass splitting effect on the 
tranverse flows is to look directly at the difference between neutron and proton 
flows 

Vr''{y.Pt)^Vriy.Pt)-V-{y,pt) 

shown in Fig. 5. 

Now the variation due to the mass splitting choice is quite evident in the whole 
Pt range of emitted nucleons. Since the statistics is much smaller at high pts 
(large error bars in the evaluation) the effect should be mostly observed in the 
Pt^^ < 0.5 range. 



3.2 Elliptic flows 



The competition between in-plane and out-of-plane emissions can be pointed 
out by looking at the elliptic flow, which has the following expression: 



V2{y,Pt) = {^^) 
Pt 



The sign of V2 indicates the azimuthal anisotropy of emission: particles can 
be preferentially emitted either in the reaction plane (V2 > 0) or out-of-plane 
(squeeze-out, V2 < 0). 

Also in this case a good check of our effective interaction choices is provided 
by some more inclusive data in the medium energy range. In Fig.6 we report a 
comparison with the FOPI — GSI data in the interesting beam energy region 
of the change of sign of the elliptic flow for protons, integrated over large 
transverse momentum contributions. The agreement is reasonable for both 
parametrizations. Also in this case we see a slightly larger overall repulsion 
(earlier zero-crossing) for the m* > m* case. 
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Our results on the pt dependence of the neutron/proton elhptic flows in Au + 
Au reactions at 250 AMeV are shown in Fig. 7. We can see that the Py 
component rapidly grows at momenta pf''' > 0.6, causing the flows to fall 
down towards negative values. At low momenta, emitted particles undergo a 
deflection in the reaction plane, more pronounced for nucleons with smaller 
effective mass. In fact, comparing the two panels of Fig. 7 we see that around 
the maximum in-plane deflection 0.4 < p[^^ < 0.8 neutron and proton flows 
are inverted. We also flnd different slopes for protons and neutrons around 
pf'^ > 0.7, depending on the relative sign of mass splitting, in agreement with 
the previous results in the transverse flows. Around ~ 1, i.e. projectile 
momentum in the cm, the elliptic flow flnally takes negative values, with slight 
differences for the two parametrizations. 

Also in this case the best observable to look at is the pt dependence of the 
difference between neutron and proton elliptic flows 

vr''{y,Pt) ^vny,Pt) -v,^{y,p,) 

shown in the Fig. 8. 

Again the difference is quite evident in the whole pt range of emitted nucleons. 
We note again that the statistics is much smaller at high pts (large error bars 
in the evaluation). Anyway the mass splitting effect appears very clearly in 
the low p[^^ < 0.5 range where the results are more reliable. 
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Fig. 6. Energy dependence of the elliptic flow in Au + Au semicentral collisions at 



mid-rapidity, |y' '| < 0.1, integrated over the ^ > 0.8 range. FOPI data with 
similar selections, see ref. [26], are given by the full triangles. Solid circles refer to 
the case m* > m*, empty circles to m* < m*. 
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Fig. 7. Elliptic flow of protons and neutrons in Au+Au reaction at 250 AMeV, 
b/b 

max — 0-5, in the rapidity interval 0.7 < |?/'''^''| ^ 0.9. 

As a general comment we stress the importance of more exclusive data, with 
a good pt selection. Indeed from Figs. 4, 5 (for Vi) and Figs. 7, 8 (for V2) 
we clearly see that more inclusive Vi{y), V2{y) data, integrated over all pts 
of the emitted nucleons at a given rapidity, will appear to be less sensitive 
to isospin — MD effects. As already discussed from general features of the 
mean field MD we expect a kind of compensation between low and high pt 
contributions. 
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p; 

Fig. 8. Pt dependence of the difference between neutron and proton transverse flows 
in Au + Au colhsions, same energy, centrality and rapidity selection as in Fig. 7. 
Solid circles refer to the case m* > m*, empty circles to m* < ni*. 
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3.3 Low rapidity results 



In order to confirm the interpretation of tlie previous flow results we liave 
performed the same analysis for a low rapidity selection < \y^^^\ < 0.25. 
Since we are now looking at particles with overall smaller momenta, the effects 
of the isospin contributions to the mean-field- MD should be reduced. 

We illustrate this point in Fig. 9, where the difference between proton and 
neutron flows are plotted as a function of transverse momentum in the low 
rapidity region. We have behaviors quite similar to the corresponding Figs. 5, 
8, but the effective mass splitting cannot be clearly observed. 

We can see some indications of a larger attraction felt by low momenta protons 
in the case m* > m*, see the p-solid line of Fig.2. This can be observed in 

the reduction of the proton transverse flow at low pf'^ , solid circles below the 
empty ones in the range pf'^ < 0.5. The effect looks interesting, although not 
much appreciable because of large uncertainty. 

3.4 Changing the stiffness of the symmetry term 

It is well known that the collective flows for asymmetric systems are sensi- 
tive to the density dependence of the symmetry term of the nuclear EOS, in 
particular at high transverse momentum [16,17,18]. 




Fig. 9. pt dependence of the difference between neutron and proton transverse 
(left panel) and elliptic flow (right panel) in Au + Au collisions at 250 AMeV, 
b/b 

max — 0.5, in the rapidity interval < ^ 0.25. Solid circles refer to the case 
m* > m*, empty circles to m* < m*. 
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mass splitting 


Zl 


Xq 




m* > m* 


-36.75 


0.11 


0.21 


m* < m* 


50 


4.449 


2. 



Table 2 

Values of the parameters zi, xq and x^ for opposite mass splitting, giving a similar 
soft simmetry term. 

For this reason it is important to be sure that signals previously discussed are 
mainly due to the n/p effective mass splitting and not strictly depending on 
the choice of the stiffness of the symmetry energy. 




Fig. 10. Total symmetry energy as a function of density for the two mass split- 
ting choices used in this work. Solid and dashed lines represent respectively the 
asy — stiff and asy — soft equations of state. 

In order to check this point we have performed some calculations using a sym- 
metry energy with a density behaviour very different from the one used before, 
but keeping the same n/p- effective mass splittings. We refer to it as asy-soft, 
see Fig. 10, compared with the asy-stijf pieyiouslj used (shown also in Fig.l 
for the potential part). Since aim is to leave the momentum dependence of the 
symmetry potential unchanged and to modify only the density dependence, 
we take the same zi values used before and choose new Xq and X3 values. In 
this way we can generate a softer symmetry energy at high density, of course 
with the same saturation value Esym{p^^) = 33 AMeV. Table 2 reports the two 
sets of new parameters. 

In Fig. 10 we show the total (potential plus kinetic) symmetry energy in the 
two cases, m* < m* (left panel) and m* > m* (right panel). Solid and dashed 
lines represent respectively the asy — stiff and asy — soft curves. 

Our results show that while the nucleon flows are depending on the stiffness 
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of the symmetry term, like in other calculations [16,17,18], the transverse mo- 
mentum behavior of the differences keeps the same sensitivity to the effective 
mass splittings. 

Fig. 11 illustrates the ra/p-flow differences as a function of transverse momen- 
tum for the same energy, impact parameter and high rapidity selections used 
before. Both transverse Vi~^{y,pt) and elliptic V2~"'{y,pt) flow differences 
show very small variations from the asy — stiff case (i.e. Figs. 5 and 8) over 
the whole pf'^ range, in particular in the region p[^^ < 0.5, of interest for the 
better statistics. 

This result represents an important indication that we are selecting experi- 
mental observables that directly probe the isospin effects on the momentum 
dependent part of the nuclear EOS. 



4 Conclusion and outlook 

We have presented some expected influence of the momentum dependence in 
the isovector channel of the effective A^A^ interaction on the heavy ion reaction 
dynamics at intermediate energies. 

For a given parametrization, dynamical effects of momentum dependence in 
the isospin channel are opposite for neutrons and protons, so that they can be 
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Fig. 11. pt dependence of the difference between neutron and proton transverse 
(left panel) and elliptic flow (right panel) in Au + Au collisions at 250 AM eV, 
b/b 

max — 0.5, in the rapidity interval 0.7 < ly^'^'^l ^ 0.9. Calculations are performed 
using the asy — soft symmetry energy. Solid circles refer to the case m* > m*, 
empty circles to m* < m*. 
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magnified by looking at the differences between neutron and proton collective 
flows. This kind of measurements requires the detection of neutron flows, but 
it is very useful to draw out some information about the relative sign of the 
effective mass splitting in a asymmetric nuclear medium. 

In spite of the low asymmetry we have found interesting effects from the study 
of Au+Au semicentral collisions at 250 AMeV. These results can provide some 
hints for future measurements, also with unstable beams. 

We suggest to look at the flows of particles of relatively large rapidity, in 
particular performing more exclusive experiments for various transverse mo- 
mentum selections. A compensation is indeed expected for observables where 
all pt contributions are summed up. This is physically understandable on gen- 
eral grounds due to the different interplay between low and high momentum 
mean field repulsion, for the various effective mass choices. 

These more exclusive flow observables, in particular the transverse Vi~^{y,pt) 
and elliptic V2~"'{y,Pt) flow differences, appear to be very selective probes of 
the momentum dependence of the isovector part of the nuclear EOS. The ob- 
served effects are not changing when we largely modify the density dependence 
of the symmetry term while keeping unchanged the momentum dependence 
of the symmetry mean field. 

Due to the difficulty of the neutron flow detection, good data with similar 
features could be obtained from the flows of light isobars, like ^H,^ He and so 
on. Exclusive data, with a pt selection, appear in any case of large importance. 

Collective flow measurements seem to represent a very nice exploration of 
the highly controversial sign of the n/p effective mass splitting in asymmetric 
matter, of large fundamental interest. Moreover related important effects can 
be expected in a wide spectrum of the nuclear physics phenomenology, from 
the structure of drip-line nuclei to the threshold production of resonances and 
particles in HIC with radioactive beams. 
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